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For stripes of hole rich hnes in doped antiferromagnets, we investigate the competition between 
anti-phase and in-phase domain waU ground state configurations. We argue that a phase transition 
must occure as a function of the electron/hole filling fraction of the domain wall. Due to transverse 
kinetic hole fluctuations, empty domain walls are always anti-phase. At arbitrary electron filling 
fraction (5) of the domain wall (and in particular for 5 1/4 as in LaNdSrCuO), it is essential to 
account also for the transverse magnetic interactions of the electrons and their mobility along the 
domain wall. 

We find that the transition from anti-phase to in-phase stripe domain wall occurs at a critical filling 
fraction 0.28 < 5c < 0.30, for any value of j < ^. We further use our model to estimate the 
spin-wave velocity in a stripe system. Finally, we relate the results of our microscopic model to 
previous Landau theory approach to stripes. 



I. INTRODUCTION 

An anti-phase domain wall in stripes is the state where 
the local antiferromagnetic (AFM) spin order parameter 
undergo a tt phase shift across a hole rich line. Such 
periodic stripe structures were experimentally found in 
doped Copper oxides and Nickel oxides [Q. Historically, 
this well accepted feature was first considered to be a nat- 
ural outcome of mean-field theory Fermi surface instabil- 
ity [^,^. Yet, no similar rigorous microscopic explanation 
was given within the frustrated phase separation picture 
Q which is currently regarded as underlying the micro- 
scopic origin of stripes in the relevant materials |^] . One 
of our goals here is to close this gap in the theory. 

In contrast with common folklore, we show that hole 
rich lines are not necessarily anti-phase domain walls of 
AFM spin domains. First, on simple general grounds, we 
argue that there must be a phase transition from anti- 
phase to in-phase domain wall configuration as a function 
of increased electron filling fraction S of the domain wall. 
We then proceed to construct microscopic t-J models of 
the local electronic dynamics which determines the re- 
sulting spin order across a hole line. 

It is frequently argued theoretically |^,^ and exempli- 
fied experimentally that the charge segregation into 
hole rich stripe lines is prior to the establishment an- 
tiphase spin domains. Therefore, the microscopic mecha- 
nism by which the hole lines enslave the spin order should 
be distinguished and considered separately. 

We develope a qualitative and quantitative microscopic 
understanding of the domain wall magnetic order. Our 
analysis accounts for the electronic dynamics both trans- 
verse and along a pre-established hole rich line between 
two AFM domains. Thus we focus solely on the mecha- 
nism which give rise to the spin antiphase domain wall 
feature, by examining the competition between anti- 
phase and in-phase configurations at given electron filling 
fraction S of the hole rich line. 

In the stripes literature it is more common to describe 
the domain wall filling in terms of the number of holes 



(below half-filling) per site along the domain wall; 
= 1- 25. 
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Experimentally anti-phase domain wall stripes in 
Nickel-Oxides have n'' « 1 (corresponding to one hole 
per site along the domain wall, or equivalently an elec- 
tron empty domain wall, 6 = 0), while anti-phase domain 
wall stripes in Copper-Oxides have n'' « i (correspond- 
ing to one hole per two sites along the domain wall, or 
electron 1/4 filled domain wall, S « 0.25). These domain 
wall filling fractions remain roughly constant over a wide 
range of dopings in the respective materials. 

We find that the transition from anti-phase to in-phase 
stripe domain wall occurs at critical filling fraction 



0.28 <6r< 0.30 



(2) 



depending on the value of ( j) . In other words, for exam- 
ple, we predict that stripe domain walls with a hole den- 
sity of one hole per three sites are always in-phase and 
not anti-phase. Appropriate numerical simulations can 
be constructed to get exact numbers beyond the limits 
of our analytical approximations. Yet, since our analysis 
indicates that the sensitivity to ( j) is rather weak, we do 
not expect the exact numerical results to deviate much 
from our predictions. 

We further apply our model to evaluate the of spin- 
wave velocity in a stripe systems, and compare with ex- 
periments. In addition, we relate the competing interac- 
tions in our microscopic model to previous Landau theory 
approach to stripes order 1^,^. Thus, we advance to- 
wards a coherent microscopic and phenomenological un- 
derstanding of the stripes structure within the frustrated 
phase-separation picture. 



A. Overview of the model construction, analysis and 
main results 

Since the paper is quite long, we here supply the reader 
with an overview of the gist of our model development 
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and main results. 

The general intuitive argument for the domain wall 
transition is the following: In one extreme case, = 1, 
where there is a hole on each site along the domain wall 
(i.e., it is empty of electrons, ^ = 0, as in Nickel-Oxides), 
it is clear that an anti-phase domain wall configuration 
is favored by transverse hole fluctuations. Now, consider 
the state of domain walls with increasing electron filling 
fraction. In the opposite extreme case where the domain 
wall is half-filled with electrons (i.e., n'^ = 0, 6 = ^), we 
should recover the undoped ordered AFM state. There- 
fore, there must be some intermediate critical electron 
filling fraction 6c of the domain wall at which there is a 
transition from an anti-phase to in-phase local AFM spin 
configuration across a hole rich line. 

The main quantitative objective of this paper is to de- 
termine the critical domain wall filling fraction Sc as a 
function of the t-J model parameters. In the process, we 
develope an advanced qualitative and quantitative un- 
derstanding of the various competing local interactions 
in the single stripe physics. Additional applications are 
discussed in sections-V and VI. 

Intuitively, increase of electron filling of the domain 
wall amounts to increase of magnetic interactions across 
the domain wall until they are strong enough to dominate 
over the charge fluctuation dynamics (due to holes). The 
possibility of a transition from topological (anti-phase) 
to non-topological (in-phase) stripes, due to increase of 
AFM interactions, was first speculated by Neto & Hone 
[0,^ (though, on not quite rigorous grounds, it was some- 
how related to the spin correlation length). 

We find that the topological/non-topological nature of 
the stripe charge wall can be completely determined by 
the local dynamics, which in turn is determined by the 
electron filling fraction of the wall (assuming fixed given 
t-J model parameters). We do not see a way by which 
interaction between domain walls, the width of the spin 
domains or any similar long length scale are significantly 
relevant. 

In section-II, we list some preliminary assumptions of 
our model of the stripe domain wall: (a) The effective 
electron dynamics is captured by a one band t-J model 
for electrons hopping between Copper sites (i.e., in which 
the Oxygen sites are integrated out), (b) It is assumed 
that the hole line, its mean position, and hole density are 
already pre-established by some higher energy processes 
(presumably phase separation and coulomb frustration). 
We hence focus solely on determining the preferred spin 
configuration across a pre-established hole line, (c) The 
spin order is determined by comparing the energetics of 
anti-phase and in-phase domain wall configurations (and 
not with the motion of dilute holes in a uniform AFM as 
was done in most previous work |p|,p^). 

In section-Ill, we introduce our microscopic model and 
start with the consideration of only transverse interac- 
tions and dynamics, (i.e., ignoring the effects of dynam- 
ics along the domain wall). The kinetic energy due to 
transverse hopping of holes favors an anti-phase domain 



wall configuration, while magnetic interaction between 
electrons favors an in-phase domain configuration (see 
figure- 1 below). These competing interactions determine 
the preferred local spin configuration. All calculations 
are done to the first significant order in {j)- Our anal- 
ysis surprisingly shows that if dynamics along the do- 
main wall is ignored then transition to an in-phase do- 
main wall would have occurred already at n'* < 2/3 
(i.e., at a density of two holes per three sites), in con- 
flict with experimental observations of anti-phase stripes 
in (LaNd) SrCuO materials with n'' « 1/2. Hence, it 
is essential to investigate the effect of kinetic dynamics 
along the domain wall, which we undertake in section-IV. 

In section-IV, we model the kinetic dynamics of the 
electrons moving along a stripe domain wall in an effec- 
tive external magnetic mean field due to its AFM envi- 
ronment. Along an in-phase domain wall there is an ef- 
fective net staggered external magnetic field, while along 
an anti-phase domain wall the net external field is zero 
on each site. We analyze two extreme limits: (a) Non- 
interacting electrons moving along the domain wall, and 
(b) Large U ^ t limit for the interaction of electrons 
along the domain wall. In both cases, the essential result 
is that kinetic fluctuations along the domain wall weaken 
the average magnetic interaction energy which favors an 
in-phase domain wall, and thus extend the stability of an 
anti-phase domain wall conflguration to higher electron 
filling fractions (i.e., lower hole densities). 

In section-V, we relate the competing interactions in 
our microscopic model to previous Landau theory ap- 
proach to stripes order ||,^ . 

In section- VI, we apply our model to evaluate the spin- 
wave velocity v±_ in a stripe state compared with vq in the 
parent AFM material (where v± is velocity perpendicular 
to the stripes). 

We argue that our results are quantitatively accurate 
well beyond the seemingly rough approximations of our 
simple models. The heart of the matter is the fact that 
our quantitative results are sensitive only to the ener- 
getic difference between an anti-phase and in-phase do- 
main wall configuration. Processes which are neglected in 
our treatment (e.g., deeper penetration of hole hopping 
into the AFM environment) have the same contribution 
in both domain wall configurations and thus only very 
weakly affect the energetic difference between them. 

II. PRELIMINARY ASSUMPTIONS 

The stripes characteristics can vary to include both 
diagonal and vertical stripes, which may be insulating 
or conducting. Therefore, the anti-phase domain wall 
mechanism should be rather simple, robust, and not too 
sensitive to the above mentioned variations. 

1. The most microscopic Hubbard type model of the 
Cu02 or Ni02 planes includes distinct Oxygen and 
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Copper (or Nickel) orbitals bands. Yet, as com- 
monly argued, we assume that the effective dy- 
namic is captured by a one-band t-J model (in 
which the oxygen sites degrees of freedom have been 
integrated out) where electrons hop directly be- 
tween Copper lattice sites. In the context of stripes 
theory, the above assumption is supported by the 
fact that correct domain wall configurations turned 
out in numerical simulations [ p2[ of one band t-J 
models. 

2. Experimental evidence and theoretical consid- 
erations jsj suggest that stripe formation is com- 
monly charge driven, i.e., that periodic hole line 
stripes form first and enslave the formation of the 
AFM spin domain. Therefore, for our purpose in 
this paper we take the hole lines to be pre-formed. 

3. Stripes were found in both Spin-i and Spin-1 doped 
antiferromagnets . A doped hole can thus corre- 
spond to a spin-0 or a spin-i site respectively. Yet, 
in both cases the hole and its dynamics are carried 
on only within the dx-i-y2 orbital band. Indeed, 
the model which we construct and analyze works 
equally well for both doped Spin-^ and Spin-1 an- 
tiferromagnets. We chose to present our analysis in 
terms of a doped Spin-i AFM (i.e., corresponding 
to stripes in a Cu02 plane). 

4. As with superexchange mechanism of antiferromag- 
netism in the undoped parent system, we argue 
that the domain wall spin order is determined by 
local interactions across the hole rich line. Hence, 
it is sufficient to consider a single stripe segment in 
isolation (see figure-1). 

5. There's a non-trivial distinction between site- 
centered and bond-centered domain walls [ p^ , 
in the sense that the spin alignment across an 
anti-phase domain wall is antiferromagnetic for 
site-centered stripe and ferromagnetic for bond- 
centered stripe. The presentation in this paper 
is conducted in terms of site-centered stripes, and 
elsewhere [|l^ we will show that the same principles 
apply for bond-centered stripes as well. 

6. Probably the main quantitative approximation in 
our model is that we treat the AFM regions be- 
tween the hole lines as if they were antiferromag- 
netically ordered. We neglect spin exchange fluctu- 
ations and the quantum nature of the AFM corre- 
lations in the rather narrow ladder geometry of the 
stripes. In other words, our quantitative results are 
rigorously valid for an Ising model approximation 
of the AFM regions. Yet, in all of our calculations 
we use a parameter e which is defined to be the 
energy difference between parallel and anti-parallel 
near-neighbor spin states (see equation (||)). Only 
in the end we substitute e = J for intuitive con- 
creteness. In principle, all the effects of fluctuations 



etc. can be incorporated into a renormalized value 
of £ without changing our results. 

7. Our quantitative results are sensitive only to the 
energetic difference between an anti-phase and in- 
phase domain wall configuration. Therefore, many 
processes which are neglected in our treatment 
(e.g., deeper penetration of hole hopping into the 
AFM environment, magnetic interaction between 
electrons on the wall) have the same contribution 
in both domain wall configurations and thus only 
very weakly affect the energetic difference between 
them. 



III. THE EFFECT OF TRANSVERSE 
INTERACTION 

Our analysis proceed in two stages; First, in this sec- 
tion, we consider solely the transverse fluctuations of a 
hole and magnetic interaction of electrons in the stripe, 
while the holes/electrons configuration along the domain 
wall is taken as static. In a second stage (section- IV), we 
consider the implications of electron mobility along the 
stripe. 

A. Model of site-centered stripe 

The average hole line position is fixed by introducing 
a chemical potential shift (/i) on particular sites [0,1). 
(In Figure-1, hole line /x-shifted sites are represented by 
dark circles). In the ground state, holes are preferentially 
situated on the ^ shifted sites. Hence, the effective local 
chemical potential shift, /i, incorporates the net effect of 
the high energy dynamics (e.g., coulomb frustrated phase 
separation [Q) which lead to the stripe formation. The 
magnitude of determines the stripe stability, and can 
be associated with the stripe creation temperature for 
the purpose of extracting an experimental estimate of 
it's magnitude. Thus, it is assumed that fi < J < t. 

An antiferromagnetic spin exchange interaction J ex- 
ists between electron spins on nearest neighbors sites. In 
order to fix the spin order, it is enough to determine the 
relative orientation of the spins immediately to the left 
and to the right of the hole line. Therefore, we include a 
transverse hoping interaction, t, only between a domain 
wall site (indicated by a dark circle in figure-la) and its 
right and left nearest neighbors in the AFM regions. As 
noted before, we argue that processes of further penetra- 
tion of a hole into the AFM regions are not significantly 
affecting our results. Hopping interaction t|| along the 
domain wall will be considered in section-IV. 

In the absence of kinetic motion along the domain wall, 
the Hamiltonian describing a single site centered domain 
wall is 
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Hi — fino_i + ^2 5^ ^-^'^ ' ^i''*' 



(3) 
(4) 



[33 ; 



co.i,. 



h.c. 



where Sj^i is the spin of an electron at column j and line 
i (where j = is the domain wall position), no.i is the 
occupation number of a domain wall site at position i 
along the wall. c\ ^ ^ and c^i ^ ^ are electron creation op- 
erators respectively to the right and left of the domain 
wall site rto,i- (t-J model projection operator which ex- 
cludes double occupancy is implicitly assumed through 
out the paper). 



B. Kinetic transverse hole fluctuations 

In the absence of hoping interaction, t, the anti-phase 
stripe groundstate (a) and the in-phase stripe ground- 
state (a') are degenerate. But, the energy of the corre- 
sponding excited states (b) and (b') differ. As a result, 
transverse hole hoping interaction removes the ground 
state degeneracy in favor of anti-phase domain wall p5[. 



Anti-phase stripe 
I- 

i o 

I • 1 I 1 
(a) (a') 



In-phase stripe 



(b) 



(b') 



FIG. 1. The domain wall sites are indicated by colored cir- 
cles. The static ground state configuration of anti-phase and 
in-phase domain walls are depicted in figures (a) and (b) re- 
spectively. The corresponding excited states (a') and (b') re- 
sult from transverse hoping dynamics of holes. Note the "bad 
bond" created by the hoping in (b'). The depicted domain 
wall groundstate segments have one electron per 5 sites on 
the domain wall (i.e., a filling fraction 5'^ = 0.1). Note the 
difference in the magnetic bonds of the domain wall electron 
between (a) and (b). 

To be general, we define e to be the energy difference 
between an antiferromagnetic "good bond" and a ferro- 
magnetic "bad bond" of two neighboring spins. Obvi- 
ously, e is proportional to J (and e = J in the case of 
Ising model of the AFM spin domains). We label by 
j^anti ^j^^ j^m bare excited states energy in the case 
of an anti-phase and in-phase domain walls as depicted 
in figure (16) and (16 ) respectively. Clearly, 



The Hamiltonian (m is thus given by the matrix 



t 



t 



t E'^ 
t E? 



(6) 



{a — anti/in) which can be diagonalized exactly, resulting 
with the ground state energy E^ 



Ef 



\/8t (l 



16*2 



(7) 



2r 

El 



for J < t 
for J > t 



The difference in kinetic energy gain due to transverse 
hoping, between an anti-phase and an in-phase domain 
wall configurations, is given by 




for J < i 
for J > i 



(8) 
(9) 



(5) 



Our result agrees with a previous large-d calculation 
(valid only in the limit J ^ t) [M. But, in the ex- 
perimental systems of interest j and thus they are 
better approximated by calculations in the limit J t 
which will be assumed for the rest of our discussion. 
From equation we conclude that, in a site-centered 
stripe geometry, the zero point transverse kinetic fluctu- 
ation of the holes favor an antiferromagnetic alignment 
of the neighboring spins. 

Note that the above result is in stark contrast with the 
conventional wisdom that hole's zero-point kinetic fiuc- 
tuations always favors ferromagnetic alignment of their 
environment p7| , [l8|| , which is based mostly on isolated 
hole models. It is a rather simple demonstration of the 
difference between collective and single hole properties in 
doped antiferromagnets. 



C. Competing magnetic interactions at electron 
fliling (5/0 

From the analysis of the above model, we conclude that 
the transverse kinetic fluctuations of holes are sufficient 
to induce an anti-phase domain wall ground state in an 
empty domain wall S = (i.e., one hole per site along 
the wall, = 1) oi stripes. 

For a domain wall with electron filling fraction of S, 26 
is the number of electrons per site along the domain wall 
(in the large U limit), and n/i = (1 — 26) is the number of 
holes in the lower Hubbard band. Each hole contributes 
a kinetic energy difference AE^™. Hence, at an arbitrary 
electron filling fraction 6 of the domain wall, the average 
transverse kinetic energy gain per site AE']^" (6) of an 
anti-phase domain wall in comparison with an in-phase 
domain wall is 
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ASf " (,5) 



kin 



(1 - 2^) 



(10) 



For each electron on the domain wall there is an energy 
difference Ai5™°^ = +e in favor of an in-phase domain 
wall, due to spin exchange interaction with the AFM en- 
vironment. At electron filling 5, the average magnetic 
energy difference (of an anti-phase domain wall in com- 
parison with an in-phase domain) per site along the do- 
main wall is 



AE' 



mag / cn _ 



\7nag 



-e{25). 



(11) 



The competition of these transverse interaction alone 
would predict a transition as a function of domain wall 
filling from anti-phase to in-phase domain wall when 



at 



AE^ {5) = AEf' + AE 



pniag 



> 



(12) 



(13) 



(corresponding to n'' < 2/3, i.e., at a density of two 
holes per three sites), in conflict with experimental ob- 
servations of anti-phase stripes in (LaNd) SrCuO mate- 
rials with 5 Ri 1/4. This conflict is resolved in the next 
section, where we account for effect of kinetic dynamics 
along the domain wall. 



IV. IMPLICATIONS OF ELECTRON DYNAMICS 
ALONG THE DOMAIN WALL 

In this section we investigate the consequences of elec- 
tron dynamics along the domain wall for the competition 
between anti-phase and in-phase domain wall configura- 
tions. As a first order approximation, we assume static 
spin configuration of the antiferromagnetic domains on 
the left and right of the domain wall. Thus, electrons 
moving along the domain wall are effectively modeled as 
a one-dimensional electron gas (IDEG) in a static ex- 
ternal magnetic field. At this mean-field level, electrons 
moving in an anti-phase domain wall experience a net 
zero external field on each site, while electrons moving 
on an in-phase domain wall experience a staggered ex- 
ternal magnetic field of magnitude proportional to the 
spin interaction strength J. Hence, the effective domain 
wall IDEG Hamiltonian has the form 
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Ho (B) = i|| ^ (4cj+i,, + h.c 
+2i?^a(-l)-'4c,, 



(14) 



-A*F^c]^Cj^ (15) 



where cr = ±1 for spin "f , J, respectively, and 



B 



for anti-phase 
for in-phase 



(16) 



When considering the competition of stripes versus 
droplets forms of local phase separation, Nayak&Wilzek 
Q proposed that a significant energy is gained by the 
increased mobility along an anti-phase domain wall. But 
we note that the kinetic energy gain in an in-phase do- 
main wall is practically the same as in anti-phase domain 
wall state. Therefore, it is essential to make a more care- 
ful analysis of the electronic dynamics along the domain 
wall in both cases. 

Below, we extract quantitative results in two limits: 
(a) For non-interacting electrons moving along the do- 
main wall, and (b) In the large U ^ t limit for the in- 
teraction of electrons along the domain wall. In both 
cases, the essential result is that kinetic fluctuations 
along the domain wall weakens the magnetic interaction 
energy gain which favors an in-phase domain wall en- 
vironment, and thus extends the stability of anti-phase 
domain wall configuration to higher electron filling frac- 
tions (i.e., lower hole densities). In particular, we con- 
clude that dynamics fluctuations along the domain wall 
allow for the establishment of anti-phase stripes with do- 
main wall filling 5 J (as in Copper-Oxides systems). 



A. Effect of motion along the domain wall for 
non-interacting electrons at arbitrary filling 

We here model the dynamics along the hole rich do- 
main wall by an effective one dimensional lattice model 
( |r^ of non- interacting electrons {U — 0). For an in- 
phase domain wall, the exchange coupling to the spins in 
the AFM environment result with an effective external 
staggered magnetic field on the electrons moving along 
the domain wall (see figure- lb). Our model of kinetic mo- 
tion along the domain wall included only band structure 
effects due to a static spin configuration of the AFM envi- 
ronment, but no dynamic scattering interactions (leading 
to finite resistivity) ]T^ . 

The anti-phase domain wall spectrum (B = 0) is 



^(anU) ^ _2tj| COS (fc„a) ~ ^f"*'^ 

27r , A^ N, 
k„a = — n ( < 71 < — ) 



(anti) rt. 

= COS 



2ti 



np 



(17) 
(18) 

(19) 



where N is the number of sites. 

For an in-phase domain wall [B ^ 0), the staggered 
field doubles the unit cell. The non-interacting Hamilto- 
nian (|l^) Hq {B ^ 0) is diagonalized by the appropriate 
Bloch states; 
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(20) 



Wj,. (fc) 



i=i 
1 



ika t 



-'^(f)Tv/(f)VcosMfc) 
cos (fc) 



The resulting energy spectrum for the in-phase domain 
waU is 




( < n < — ) 



2nF 



(22) 



(23) 



In this context we comment that for a Hubbard model in 
an external staggered magnetic field, unlike the case of 
staggered charge potential/interaction, there is no charge 
gap opening at 1/4 filling For simplicity, the rest 

of the formulas are written for the case where the odd- 
sites are with the magnetic field anti-parallel to the spin. 
Moreover, we remind the reader that we assume B ~ 
J < Weh below half filling, (i.e., i|| cos (fc) > 0), 



-B + y^W+tfcoPlk) 



1 - 



i|| cos (fc) 
B 



til cos (fc) 



B 



til cos (fc) 



(24) 



The gain in magnetic energy per fc-state is due to the 
difference between the occupation probability of odd and 
even sites for a given electron spin, (this is the main 
difference between IDEG in an in-phase vs. anti-phase 
domain wall), 



'^2j,(7^2j,a- 



C2j-l,a'^2j-l,a j ^ 



(25) 



^ 2 



1 



B 



B 



fj| cos (fc) 



t|| cos (fc) 
B 

i|l cos (fc) 



Therefore, for filling up to kp^ the magnetic interaction 
energy gain for an in-phase domain wall in comparison 
with an anti-phase domain wall is given approximately 
by, 



^^mag _ ^B 



7t5 



B 



-ttS h COs(fc) 



dk 



(26) 



= 2B 

_ J 

4 



In 



In 



1 -I- sin (ttS) ' 
1 — sin {tt6) ^ 
sin (ttS) ' 



1 — sin (ttJ) 



(27) 
(28) 



(Note that the limit t|| of the previous section cannot 
be recovered since we used approximations based on t|| 3> 

We are now in a position to give a first analytical esti- 
mate of the critical transition point filling 6c determined 
given by (using (^ and (^)) 



= AEi_ {S, 
J 



(1 - 2Sa) + 



J 

4 



J 



In 



1 + sin {nSc) 
1 — sin [ttSc) 



(29) 
(30) 



The geometrical solution for j — jq is plotted 

in figure-2 




FIG. 2. geometrical solution for j — Jq 

The main qualitative result is that electron hoping fluc- 
tuations along a partially filled domain wall inhibit the 
effective magnetic interaction energy across the domain 
wall. As a result, an anti-phase domain wall ground- 
state configuration can be obtained beyond S = ^, de- 
pending on the value of J/t. For ^ < j < jq we find 
0.3 < 5c < 0.4. Note that the critical filling fraction Sc is 
not very sensitive to the value of J/ 1. Most importantly, 
it is above the value 6 « 0.25 of anti-phase stripes in 
(LaNd) SrCuO systems M. 



B. Large U ^ t model, and a second estimate of the 
critical domain wall filling 5c 

In this subsection we examine the competition between 
anti-phase and in-phase domain wall configurations for 
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strongly interacting electrons, U ^ t, along the domain 
wall at arbitrary filling fraction. 

The exact ground-state energy of the Hubbard model 
(|lj) is well known for the case B = from Bethe ansatz 
solution. Unfortunately, we do not know of an estab- 
lished solution for the groundstate energy of a one di- 
mensional Hubbard model in a staggered magnetic field 
B ^ 0. In particular, when B = the Hubbard inter- 
action in momentum space takes the form [/nfein^i, and 
the U = oo limit is effectively captured by imposing a 
restriction of one electron per /c— state. But such a sim- 
ple representation does not exist in terms of the Wj^a {k) 
basis (|2^) in a staggered external magnetic field. There- 
fore, we here develope a way to approximate the two 
competing wall configuration energies at the same level 
of approximation. 

We wish to tress that the qualitative effect of longitudi- 
nal hoping fiuctuations - to inhibit the effective magnetic 
interaction energy across the domain wall - was already 
established in the preceding subsection. The purpose of 
the current subsection is to investigate the quantitative 
effect of including large Hubbard interactions between 
the electrons. Though we indeed resort to approxima- 
tions at several stages of our model and calculations, we 
do capture the substantial effect. Once the core physics 
established, the degree to which our final numerical val- 
ues are exact can be checked (and improved) in numerical 
simulations on finite systems. 

Our main result is that the large Hubbard interaction, 
U, delimitate the critical filling fraction to be in the nar- 
row range 



0.28 <Sr< 0.30 



(31) 



for any value of j < ^ . 

Below, we explain our modeling approach and calcula- 
tion. We first apply the model to examine the simplest 
case of 1/4 filled domain wall, and show that always an 
anti-phase is obtained. Then, we apply the model to re- 
estimate the critical filling fraction for a transition to an 
in-phase domain wall, and thus establish the result noted 
above in equation (pn|). 



1. Model of characteristic unit cells 

The unit cell of a stripe domain wall comprise of two 
lines as depicted in figure-3. Our modeling idea in this 
subsection is to solve exactly the Hamiltonian of iso- 
lated prototypical unit cells, and then approximate the 
full stripe as assembled of a collection of such unit cells. 
Such an approximation is building on our findings in the 
previous subsection, where we argued that it is not the 
global longitudinal mobility, but only the local kinetic 
fiuctuations (in the occupation of odd/even numbered 
sites along the wall) which distinguish anti-phase and in- 
phase domain walls. 



As shown in figure-3, in the L'^ = oo limit, there are 
only three possible electron occupation numbers of a unit 
cell: one, two or zero electrons on the domain wall sites. 
In our model approximation, we ignore the hopping from 
one unit cell to another. As we shall later explain, for 
filling fractions 5 > j it will suffice to consider only unit 
cells with one and two electrons on the domain wall sites 
(i.e., as in figure-3(al,a2 and bl,b2)). 



Anti-phase stripe 

al' 




bl 



b2 



b3 



In-phase stripe 

bl' 



tI 1 bl' 11 j 
111 IM 



FIG. 3. Characteristic unit cells of a stripe domain wall. 
Any domain wall groundstate can be constructed by prop- 
erly assembling these unit cells. For domain wall groundstate 
with filling fraction 5* > 1/4, unit cells (a3) and (b3) can be 
neglected. 

We are interested only in the effect of hopping dynam- 
ics t|| of electron along the domain wall. Therefore we 
examine first the unit cells with one hole and one electron 
on the domain wall sites. There's spin exchange interac- 
tions of strength J between nearest-neighbor electrons. 
Yet, to zeroth level of approximation, the spin config- 
uration is assumed to be static (apart from t|| electron 
hoping). Thus, we consider only the electron configu- 
rations of an anti-phase (al,al') and in-phase (bl,bl') 
domain wall unit cells. 

Consider the case of unit cells with a single hole and 
spin- J, electron hopping between the domain wall sites. 
We observe that while the anti-phase domain wall states 
(al,al') are degenerate, the in-phase domain wall states 
(bl and bl') have an energy difference of magnitude 2e 
(where e ~ J is the energy difference between paral- 
lel and anti-parallel neighboring spin states). Thus, the 
Hamiltonian of a single electron in an anti-phase domain 
wall unit cell is 



Tranti 
holc-c 



-t 



-tu 







(32) 



with groundstate energy —t, while for an in-phase domain 
wall 



1 



with 



holc-coll \ 

groundstate 



£2+^2 



tn + 



2 (t||) 



(33) 
energy 

(Note that the mag- 



netic interaction energy is already fully accounted for in 
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the eigenvalues, so it shouldn't be counted again). Thus 
we find that, for a stripe unit cell with one hole and one 
electron on the domain wall sites, the ground state mag- 
netic interaction energy difference between an anti-phase 
and in-phase domain wall is 



A rpmag 



T^anti 
-'^holo-ccll 



hole-cell 



(34) 



= -ill 



2 (*li) 



II 

for t|| 
for ill 



> J 
^ 



The above needs to be divided by 2 in order to get the 
energy difference per site, since the unit cell has two sites 
along the domain wall. In the static limit t\\ 0, we re- 



cover equation ( |ll| ) of section-Ill, ^ ' 2' 

But for the rest of the paper, we focus on the experimen- 
tally more relevant limit t|| ^ J. 

We now proceed to demonstrate first the essence of our 
approach for the simplest case of 5 = ^ filling. 



2. The special case of 1/4 filling 

In section-Ill, we have shown that dynamics of only 
transverse interactions (hole hopping and Heisenberg in- 
teraction) will not suffice to stabilize an anti-phase do- 
main wall for S > 1/6. Therefore, by investigating the 
case of (5 = 1/4 we can already answer the question of 
principle about the effect of electron dynamics along a 
stripe domain wall. Moreover, it is a case of particular 
interest for stripes in HTc, which are found with a filling 
fraction very near 5 = 1/4. 

In the limit of large on-site interaction U ^ t, at 
(5 = 1/4 there is effectively one electron for every two sites 
along the wall. Moreover, to mimic the expected effect 
of coulomb interactions, it makes sense to supplement 
the model with a repulsive interaction also between near- 
neighbor holes (i.e., a, t — U — V model along the domain 
wall) pO| . Thus, in the absence of kinetic hoping dynam- 
ics along the domain wall, the groundstate of the one di- 
mensional electron gas is a periodic charge density wave 
(CDW) with one electron per every second site. Even 
with the inclusion of kinetic fluctuations, the same CDW 
still dominate the electron correlations at low tempera- 
tures. As a result, we argue that the occurrence of unit 
cells with an occupation of zero or two electrons on the 
domain wall sites (as in figure-3(a2,a3)) are rare events 
and thus their contribution may be safely neglected. 

Therefore, a two-line stripes segment with one hole and 
one electron on the domain wall sites can be regarded as 
the characteristic unit cell along a 1/4 filled domain wall 
(as in figure-3(al,bl)). By solving such unit cell model 
exactly, we are able to account for the effect of local 
kinetic fluctuations along the domain wall. Moreover, 
since we are interested only in the difference between 



anti-phase and in-phase wall, and since the net magnetic 
field averages to zero along both anti-phase and in-phase 
wall, it is intuitively expected that the significant dif- 
ference is not in the net electron mobility but mostly in 
the local kinetic fluctuations (Though magnetic scatter- 
ing from AFM environment, which we neglect, may also 
differ). 

Let us recall first the result of section-Ill. In the ab- 
sence of hopping along the wall, the average transverse ki- 
netic energy per site of an anti-phase domain wall in com- 
parison with an in-phase domain wall was (using equation 
®) 



AE 



kin 



s = 



1 



kin 



N 



and the competing transverse magnetic energy due to 
Heisenberg interaction with the AFM environment was 
(using equation (11)) 



TV 



= +e (26) = 



Hence, it was the in-phase domain wall configuration 
which had the lowest energy ground state. Now, with 
the inclusion of the effect of the kinetic fluctuations along 
the domain wall, the average magnetic interaction energy 
is modified and (using ( p4[ ) instead of (^ij)) the average 
energy difference per site is 



— = AE^'' 
N ^ 



2 ^-^holc-c 



AEjS = 1/4) 
N 



e \ e 

Tn 



< 



(35) 
(36) 



where N is the number of sites along the domain wall, 
(The factor 1/2 in (|3^) is due to having two domain wall 
sites per unit cell). Therefore, we find that due to the 
kinetic fluctuations along the domain wall, the magnetic 
interaction energy gain of an in-phase domain wall con- 
figuration is reduced enough to make the anti-phase do- 
main wall more favorable at 1/4 filling for any value of 
J ~ e < til . 



3. Transition at domain wall filling Sc 

We now attempt to determine the contribution of the 
effect described above at arbitrary electron filling 2(5 > i , 
and thus estimate the transition point. 

As argued before, the IDEG correlations on the do- 
main wall are dominated by a CDW correlations where, 
at electron filling 2S > 1/2, the occurrence of two near- 
neighbor holes is statistically negligible compared with 
other configurations in the ground state. Therefore, the 
two-sites unit cells along the domain wall are basically 
only of two kinds: (a) There are (1 — 26) unit cells con- 
taining one hole, which are the type analyzed in the pre- 
vious subsection, (b) There are ^ — (1 — 26) = (26 — ^) 
unit cells containing no holes. 
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To estimate the average energy contribution per site, 
we can use the resuh of the previous subsection for the 
hole cehs, only with the added factor (1 — 26); 



1 



-ptnag 



(1 - 25) 



e 








1 - 




4 





(37) 
(38) 



and add the contribution of magnetic energy difference 
per site due to unit cells with no holes, 



1 



AE" (S) = -E 



1 



1 



(39) 



The preferred groundstate configuration is determined by 
AE {5) 



N 



AE" (5) + AE" {S) 



(40) 



1 



(1-26)+ [25- -]e. (41) 



A transition from anti-phase to in-phase domain wall oc- 
curs when 

AE (5) > 

Therefore, the critical filling fraction 5c is given by 



'—2 

3 _ 
2 



(1 - 25c) + 2[25c-^]e (42) 



(43) 



The predicted 5c for several values of | is given below. 



£ _ 1 
t " 3 



= 0.285 71 



I I = ^ 1 = 0. 293 1 



<0.3 



We find that 0.28 < 5c < 0.3, (for j < and is only 
weakly sensitive to the value of j . 



V. CONNECTION WITH A LANDAU THEORY 
OF STRIPES ORDER 

The properties of a general Ginzburg-Landau free en- 
ergy (|j) of stripes were previously investigated in [^|| . 
Specifically, an ordered stripe phase is a unidirectional 
density wave which consists of coupled spin-density wave 
(SDW) and charge-density wave (CDW) order parame- 
ters. 



J^q,k = ris |S^| -I- r2s 

+7 [Sts.^,rp;~ + c.c 





2 


q+k 





(44) 



where is the complex-valued CDW amplitude with 
wave vector fc, pi = and similarly S^- is a complex 
vector amplitude of the SDW. The quartic (and higher 
order) terms required for stability are omitted. 

Zachar et al. have considered the phase transition 
between a stripe phase and a high-temperature disor- 
dered state, as involving only a single SDW wave vector 
(i.e., SI = S-_|_gj. The existence of the cubic 7 term, 
coupling these two order parameters in the Landau free 
energy, drives the period of the SDW to be twice that of 
the CDW, and the absence of any net AFM ordering is 
equivalent to the statement that the stripes are topolog- 
ical. 

In contrast, it was shown by Pryadko ct al. that 
the same sort of cubic term in a Landau theory of the 
transition from a homogeneous ordered antiferromag- 
netic phase to a stripe ordered phase produces a state 
in which the antiferromagnetic magnetization does not 
change its sign between the domains, i.e. the stripes are 
non-topological . 

When investigating the transition from a well- 
developed antiferromagnet with a modulation vector 7? = 
(tt, tt) to an incommensurate modulated phase, we must 
account for both the original AFM order parameter 
(which cannot be assumed small), and the spin density 
wave , r. The most relevant terms in the Landau ex- 
pansion of the effective free energy are then 



T = ris |Sf?| 







^ 1 |2 


'"2s 


^a+k 





(45) 



where a finite ris < is assumed as given. This expres- 
sion implies that an instability in either spin or charge 
sector generates both spin- and charge-density waves at 
the wave vectors q — ir + k and k, respectively. Near the 
transition the magnitude of the incommensurate peak is 
necessarily much smaller then the commensurate AFM 
modulation |S-^^:| ^ |Sff|. It is easy to see that this 
corresponds to in-phase domain walls; The derived rela- 
tionship between q and k implies that the periods of spin 
and charge modulation are equal. 

Our analysis in this paper supply microscopic insight 
to the Landau theory results. First, we find that the 
possibility of both topological and non-topological stripe 
phases can indeed be realized within a single microscopic 
model (as anti-phase and in-phase domain wall ground- 
states respectively). 

Second, we find that non-topological stripes (corre- 
sponding to in-phase domain wall state) are indeed es- 
tablished due to enhanced antiferromagnetic interactions 
(as first speculated by Castro-Neto Q). Yet, thus far the 
connection is more heuristic than rigorous. 
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Furthermore, after analyzing the effect higher order 
derivative terms in a gradient expansion of the Ginzburg- 
Landau free energy, Pryadko et al. argue that §: When 
there is no frustration, topological stripes are not estab- 
lished in the ground state. However they speculate that 
frustration, such as competing first and second neigh- 
bor interactions, or opposite sign terms in the gradient 
expansion of the Ginzburg-Landau model {i.e. below a 
Lifshitz point), can stabilize topological collinear domain 
walls. In other words, topological stripes are a conse- 
quence of physics on short length scale, and they do not 
appear in a theory that considers only long-distance or 
low- energy physics. 

Pryadko et al. speculate that the missing frustrating 
ingredient may be due to inter-stripe interactions such 
as long range coulomb interactions which are expected 
to exist between the charged domain walls. Our analysis 
in this paper demonstrate that such interactions are not 
required. We have shown that local single stripe dynamic 
interactions, (albeit involving finite strength competing 
interactions, i.e., expectantly beyond the reach of pertur- 
bative gradient expansion) are sufficient to realize both 
alternative possible topological magnetic order across the 
charge line. 



VI. ESTIMATE OF SPIN- WAVE VELOCITY 

In this section, we use our microscopic domain wall 
model to estimate spin-wave velocity reduction 



v± 



(46) 



where v±_ is the spin-wave velocity perpendicular to the 
stripes (i.e., along the modulation direction), and Wq is 
the spin-wave velocity in the undoped parent antiferro- 
magnetic material. Equivalently, it is the predicted spin- 
wave velocity anisotropy in stripes (where is the same 
as the velocity along the stripes) . Castro-Neto & Hone |^] 
were the first to point out that within the stripes model, 
one of the main effects of the hole rich lines would be to 
weaken the effective exchange interaction between spins 
on either side across the domain wall. Thus, a cardi- 
nal parameter in the model is the effective weakening 
factor of the magnetic interaction Jw across the stripes 
domain wall as compared with the Heisenberg interaction 
J within the hole free AFM regions (and hence along the 
stripes). 



J 



(47) 



The key new ingredient is our analytical estimate ofaj. 
Consequently, while previous model estimates always in- 
volved fitting of free parameters, we are able to estimate 
analytically the spin-wave velocity anisotropy (without 
any free parameters). 




Neto&Hone [Q proposed an effective anisotropic 
Heisenberg model for describing the low temperature 
spin dynamics in a striped Cu-0 plane. Within the 
anisotropic Heisenberg model 0, they derived the spin- 
wave velocity anisotropy 



(48) 



i.e., ay ~ V2^- fitting other consequences of their 
model to experimental results, Neto&Hone extracted the 
value 

aj « 0.01. 

The fitted value of aj, together with ( ^8|) hv± « 
50meV—A. As discussed by Tranquada et al. |^l[], this 
value is much too small to be compatible with inelastic 
neutron scattering experiments. It is further incompati- 
ble with the range of energies over which stripes incom- 
mensurability is observed. 

Aiming specifically to model the interaction between 
narrow stripe AFM regions, Tworzydlo et al. made 
predictions based on models of coupled narrow spin lad- 
ders (e.g., 2-leg or 3-leg ladders). The resulting expres- 
sions for the spin wave velocity anisotropy were 



ay 



2a J 
1+Qj 



for 2-leg ladder 
, , „ for 3-leg ladder 



3a J 



(49) 



In terms of 2-leg ladder model, Tranquada et al. pro- 
claimed to fit the experiments well with a value of 



givmg 



aj « 0.35. 



0.72 



In principle, the ladder model results are quite sensitive 
to the width of the stripes (which affects the ladder spin 
gap magnitude). Yet accidentally, the above value of 
aj « 0.35 would give a^ ~ 0. 78 for 3-leg ladder, i.e., 
practically the same as for 2-leg. 

Our microscopic model is based solely on the local 
physics in the neighborhood of the hole line. As such, 
it is insensitive to details of the width of the AFM spin 
domains. We argue that, by definition, the energy differ- 
ence between alternative spin configurations across the 
domain wall is a measure of the effective spin interaction 
across the domain wall in any appropriate low energy the- 
ory. Therefore, using our microscopic model calculation, 
we suggest that 



Jw = -^E (d) ^ - 
were AE (S) is given in equation (^) 



(50) 



AE{S) 



(£)](l-2.) + 2(2.- U 
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(where as discussed previously, e « J). Thus we can 
calculate aj and then, following the anisotropic Heisen- 
berg model reasoning of Castro-Neto & Hone, use (^) 
to calculate the spin-wave velocity anisotropy. 

For the particular case of stripes observed in LaNdSr- 
CuO § 5 w 1/4, we obtain 



(51) 





_ J 










1 - 


m 




4 







Jw ~ A-B 

Substituting ^ > j > jq: we obtain 

0.17 < aj{SK 1/4) < 0.23 



(52) 



Importantly, note that aj (6 sa 1/4) is rather insensitive 
to the variation of j < i. (where 6 ^ 1/4 is the sup- 
posed electron filling fraction of the stripe domain walls 
in Copper-Oxides). 

Using equation (Esl) thus give an estimate of 



0.29 < a„ < 0. 34 

or using equation ( ^9| ) give an estimate of 

0.67 < Oft, < 0.71 for 2-leg ladder domains 
0.74 < < 0.78 for 3-leg ladder domains 



(53) 



(54) 



which fit remarkably well with the experimentally de- 
duced spin wave velocity inhibition 0.60 < a^, < 0.72 in 
doped Copper-Oxides compared with the undoped par- 
ent AFM material |^^. Note that our estimate is ob- 
tained by directly substituting the experimental value of 
vq into the analytical results for aj (using our model 
equation (^)) and for (from the ladder models ||2|] 
or the anisotropic Heisenberg model Q]), without any fit- 
ting of free parameters. 
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